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Abstract. Corrected versions of the numerically invariant expressions for the 
affine and Euclidean signature of a planar curve given in [|lj are presented. The 
new formulas are valid for fine but otherwise arbitrary partitions of the curve. 
We also give numerically invariant expressions for the four differential invari- 
ants parametrizing the three dimensional version of the Euclidean signature 
curve, namely the curvature, the torsion and their derivatives with respect to 
arc length. 



1. Introduction 

The concept of signature was introduced by Calabi et al. M as "a new paradigm 
for the invariant recognition of visual objects." For the case of planar curves, the 
signature is defined as follows: if a(s) is a smooth curve in M 2 parametrized by arc 
length s and if G is any finite dimensional transformation group acting transitively 
on K 2 , then the signature curve S of a(s) with respect to G is given parametrically 
as (k(s), k s (s)), where k is the G-invariant curvature and k s its derivative with 
respect to arc length. One of the consequences of a theorem proved by Cartan [|| 
is that S fully determines the given curve a modulo G, provided that k is never 
zero. 

The signature curve can therefore be used to program a computer for recognizing 
curves modulo certain group transformations. However differential invariants are 
high order derivatives and hence very sensitive to round off errors and noise. The 
idea of writing a numerical scheme in terms of joint invariants was also introduced 
by Calabi et al. Q. Hoping to obtain less sensitive approximations, they suggested 
to find numerical expressions for k and k s in terms of joint invariants. A joint 
invariant of the action of a group on a manifold is a real valued function J which 
depends on a finite number of points x±,...,x n of the manifold and which remains 
unchanged under the simultaneous action of the group G on the point configuration, 
i.e. J{x\, ■■■,x n ) = J(g-x\ , g-x n ) V<? € G. For example, the Euclidean distance is 
a joint invariant of the action of the Euclidean group on R 2 . Expressing differential 
invariants in terms of joints invariants results in a G-invariant numerical approx- 
imation. In their paper, Calabi et al. proposed numerically invariant expressions 
for k and k s for two specific group actions, namely the proper Euclidean group and 
the equi-afhne group. 

But contrary to their claim, the expressions given for k s are not convergent for 
arbitrary partitions of the curve. In the next section, we give correct formulas for 
approximating k s and explain why the old ones do not work in generaL. In order 
to prove our claims, we also compute the resulting numerical signatures in a prac- 
tical example and with different partitions. We then compare with the signatures 
obtained with the old formulas. 
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Cartan's theorem also provides us with a way to characterize curves in IR 3 mod- 
ulo group transformation. The generalization of the signature curve is a curve in 
M 4 determined by four differential invariants: the G-invariant curvature k, its de- 
rivative k s with respect to G-invariant arc length s, the G-invariant torsion r and 
its derivative t s with respect to s. For practical applications, we are interested 
in the case where G is the proper Euclidean group. Following the example of p|, 
we have found approximations for k, k s , t and r s in terms of the simplest joint 
invariant of the action of the Euclidean group on R 3 : the Euclidean distance. The 
results are given in section 3. That section also contains the results of numerical 
tests performed on a space curve with different parametrizations. 



2. Corrections for the Case of a Planar Curve 

In principle, one must keep track of the order of the approximation when ma- 
nipulating an approximation. This is especially true when trying to approximate a 
derivative using an approximate expression. For example, if c(t) = c(t) + 0(t) is a 
first order approximation for c(t), then for < h,t <C 1 

5(t) - 5(0) = c(t) + Q(t) - c(0) ^ c/(()) + Q{tQ) 

is in general not an approximation for c'(0). 

By keeping track of the order of approximation throughout the computations 
performed in jy, one finds out that the expressions given for n s will converge only 
for a 2 nd order approximation of n. In the case of a regular partition of the curve, 
the approximations used for the curvature are of second order, as one can see 
from equation (|l]) below. This fortunate fact explains the correct numerical results 
presented in Jy . However for a generic partition there is no guarantee of success 
and failure is likely as we show by example. 

In the following we present the corrected formulas together with their justifica- 
tion for the Euclidean and affine group action.. 



2.1. Euclidean Group Action. Let P;_2, Pj-i, Pi, Pi+i, and Pj+2 be 5 consec- 
utive points on a smooth planar curve. Denote by d(Pj,Pk) the Euclidean distance 
between the points Pj and Pfc. As illustrated in figure 1, let 

a:=d(Pi_i,Pi), c:=d(P_i,P+i), e := d{P u P i+2 ), g := d(P_ 2 , P-i) 
b:=d(P h P i+1 ), d:=d(P i+1 ,P i+2 ) f :=d(Pi_i,P i+2 ) 
It was shown in [GJ that k{P; t ) — where A denotes the area of the triangle 

with sides of length a, b and c, is a good approximation for the Euclidean curvature 
at Pi . In fact, the following expansion has been proven to be valid for small a and 
b. 

(1) «(Pi) = K(Pi) + ^(b-a) Ks (Pi) + 0(2). 

Observe that k depends only on the distances between the three points. Since the 
distance is an Euclidean invariant, we thus have a Euclidean numerically invariant 
expression for the curvature. 

When a ^ b, the approximation obtained is only of first order. Therefore care 
must be taken when approximating the derivative n s (Pi). The formulas used in [pi 
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are 



(2) 

(3) a s ,i{ p i) 



b 

K(P i+1 ) - K(Pi-l) 



c 

We would like to have expressions that do not assume a priori conditions on the 
parametrizations. Hence we propose the following numerically Euclidean invariant 
expression for n s {B): 



(4) S.,3(Pi) := 3- 

(5) K s A p i) 



HPi+i) - 




a + b-\ 


-d ' 


K(Pi+l) - 


■ Z(Pi) 


a + b- 


f d 


K(Pi+l) - 


K(Pi-l) 


2a + 2b - 


^d + g 



2 a+b+d 2 a+b+g 

(6) K s , 5 (Pi) := 3- 

Claim 2.1. Formulas (^j, (Q) and (j^J converge to K s (Pi) as a, b, d, g — > 
Proof. By ([j]) we have 



and 



S(P<) = K(Pi) + i(6 - o)k.(P<) + 0(2) 

«(P 1+1 ) = k(P 1+1 ) + ~(d - 6)k,(^+i) + 0(2). 

Assume that a=a(s) is a smooth curve parametrized by arc length, that P s ; = a(0) 
and that Pi+i = a(h). We can write 

k(P(+i) = k(P) + /iK 5 (P) + 0(2) « K (P) + bK s (Pi) + 0(2) 
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and 



Thus 



K s (P i+1 ) = K s (Pi) + 0(1). 

h(P i+1 ) - k(Pi) = ^{d - 2b + a)n s (Pi) + b Ks {Pi) + 0(2). 
Rearranging, we obtain 

z- m i; + T i) =^)+o(i) 

a + b + a 

Therefore (fij) is valid. 

Equation (|J) is just a symmetric formula obtained by averaging (f|) with its 
reflected version. We obtain @ by expanding k(Pi + i) — £(Pj_i) about Pi in a 
similar way. □ 

Figures 2-8 show that our formulas are just as good as (||) and (|^) when the 
parametrization is very regular. But when the parametrization is not, we see that 
(||) and (||) diverge from the exact solution, while (^), (|J) and (0) remain closer 
to it. Thus we believe that our formulas are better for a generic partition of a. 
Observe that formula suffers from a small bias due to its asymmetry. This 
phenomenon was also observed for formula (^|) in Q . 

2.2. Affine Group Action. Let a be a convex smooth curve. Let Pq, Pi, P2, 

P3, Pi be five consecutive points on a. Denote by [ijkl] the signed area of the 
parallelogram whose sides are Pi — Pj and Pk — Pi and by [ijk] the signed area of 
the parallelogram whose sides are Pi — Pj and Pi — Pk . Define T and S at P2 by 



4-T(P 2 )= 11 [Imn], 

0<Km<n<4 

A-S(P 2 ) = [013] 2 [024] 2 [1234] 2 + [012] 2 [034] 2 [1324] 2 

- 2[012] [034] [013] [024] ([123] [234]+ [124] [134]). 

Observe that the expressions for T and S are equi-afhne invariant. The following 
is a numerically invariant approximation for the affine curvature at P^. 

o "( p -> - 5^ 

Indeed the following expansion has been shown to be valid for sufficiently close 
points 0] 



(8) ki 

where Lj denotes the signed affine arc length of the conic from Pi to Pj and the 
higher order terms are quadratic in the Ly's. Again, since this is a first order 
approximation, care must be taken when approximating the derivative of k with 
respect to affine arc length s. 
The formula used in [fil is 
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(9) K s = = , 

Si + Si-i 

where Sj is the triangular approximation to the afhne arc length from Pj to Pj+i- 
This formula would be valid only if k were a second order approximation of k. How- 
ever, this is not the case unless Lj 1S small compared to all of the L/s. We 
would like a more general expression. Hence we propose the following numerically 
affme invariant expression for n s which does not assume a priori conditions on the 
paramctrizations : 

(10) «i = 5 .„ ^k^H _ 

Si— a + 2Si_2 + 2Si-i + 2Si + 2Si±i + S,_ % 

The convergence can be proved by an argument similar to the one used in the 
Euclidean case. 

Figures 9-12 show the results obtained with our formula using different partitions 
of a curve. Comparison is made with the signature curve obtained with (^). Just 
like in the Euclidean case, we see that equations (||) and ( |l0| ) seem just as good 
when the partition of the curve is regular. But when the partition is not regular, 
(||) diverges from the exact solution while (|o|) remains valid. 

3. Numerically Invariant Euclidean Signature for a Space Curve 

The signature curve 5 of a space curve is given by S — (k, k s ,t,t s ). In the 
following section, we present some numerically Euclidean invariant approximations 
for 5 for the case of the Euclidean group action on R 3 . This is achieved by expressing 
each of the differential invariants constituting S in terms of Euclidean distances, 
which is the simplest joint invariant of the Euclidean group action on R 3 . 

3.1. Numerically Invariant Expressions for n and k s . Let Pi-i, Pi and Pi + 1 
be three consecutive points on a space curve and let a, b and c be their mutual 
distances as illustrated in figure 1. Denote by A the area of the triangle with sides 
a, b and c. The curvature of a space curve is defined exactly the same way as the 
curvature of a planar curve: it is the inverse of the radius of the oscullating circle 
at p. So we expect that k = ±^ is a good approximation for k at i^. Indeed, a 
Mathematica routine using the local canonical form of a curve in R 3 gave us the 
following Taylor series expansion. 

{b ft) , / 2 | l2\/^ ss kt^ ^ uf^ s 



k s + (a 2 + b 2 )<^- - —)- a6(— + — ) 
3 v n 12 36 ' v 12 36 ' 

+ 6 3 -a 3 )( ! + _!££.) 

v A 40 570 45 60 ; 

2 2 

+ (a 2 b-ab 2 ){ — t + + — J. + _££!) + ... 



60 180 180 60 
Note that if we set r = 0, then our result agrees with the 2-dimensional case. Also 
observe that r does not appear in the first two terms of the expansion. Therefore the 
expressions for k s given in the planar case ((|j), (0) and (Bh) are valid numerically 
invariant expressions for k s in the three dimensional case. 
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3.2. Numerically Invariant Expressions for r and r s . The Euclidean invari- 
ant torsion r is defined as the derivative of the angle of the oscullating plane with 
respect to arc-length. Let — 1 <C — 8 < Q < ei < e% <C 1. Suppose P,_i = a(— 6), 
Pi = a(0), Pj_|_i = a(ei) and P;+ 2 = afa) are four consecutive points on a space 
curve a. Let a, 6, c, d, e and / be their mutual distances as illustrated in figure 
1. Denote by H the height of the tetrahedron with sides a, b, c, d, e and / with 
respect to Pj+2 and write Ai x i 2 i 3 for the area of a triangle with sides of length l\, 
li and ^3. 

We propose the two following numerically invariant expressions for the torsion 
at P: 



(12) t - 3 Hb 

2 fA e bd 

Observe that ([ll]) and ( |l2"| ) both use the minimal number of points (i.e. four) for 
approximating r. One might argue that this induces an asymmetry in the numerical 
results. However, this is something easy to fix in a practical situation (for example 
by averaging the results obtained in the two directions). 

We justify (|l"l| ) as follows: ft is known that t = ""'""j-'"- ; In order words, 
|t| is the component of a sss which is perpendicular to the oscullating plane. We 
start by writing down a finite difference scheme for a sss (for justifications, see [|j 
§5.4). 



(13) 



a(t2) — Oi(ti) a(ei) — a(0) a(ei) — a(0) a(0) — a(— 5) 



,(())« 6- 



e 2 - ei 




ei 


ei 




5 




£2 











£2 



When ei and i5 approach zero, the plane defined by Pj_i, P, and Pj+i approaches 
the oscullating plane. More precisely, if 5 and e\ are small enough, we find (from the 
local canonical form of the curve) that the equation of the plane passing through 
a(—S), a(0) and a(e) can be written as 

Ax + By + z = 

with A = ^ei<5 + 0(3) and B = -f (ei - 5) + 0(2), Here, x, y and z are the 
Frenet frame coordinates. The normal vector to this plane is ,[ A ,' B j}) which is 
a first order approximation of (0,0, 1), the normal vector to the oscullating plane. 
We have 



kt =a sss ■ (0,0, f) 

(A,B,1) 
asss VA 2 + P 2 + f 



•0(1). 



But 



a(-5) ■ {A, B, 1) = a(0) • (A, B, 1) = a(e x ) • (A, B, 1) = 
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and 



Therefore 



(A,B,1) 
V ' VA 2 + B 2 + 1 



±6H 



(e 2 - ei)e 2 {5 + e 2 ) 
6H 



W± def 

From the previous section we know that k rs ±k = 4A ° bc and therefore k w ± 

There is an interesting similarity between (|llj) and k. This is easily seen if we 
write k in a slightly different way. By basic geometry argument we can say that 
k = ±2| where h is the height of the triangle with sides a, b and c. Thus kr can 
be looked at as a three dimensional version of k. This is one reason why we like 
( p~T| ) . Another reason is that it provides us with an easy visual understanding of the 
torsion similar to the understanding we already have of the curvature. It is indeed 
easier to evaluate the height and sides of a tetrahedron together with the radius of 
the oscullating circle than a derivative. It makes the evaluation of the torsion from 
a picture of the curve very intuitive. 

For finding (|l2|), we started from the following definition: if P = a(t) and 
Q = a(0), then r = limp^g |, where is the angle between the oscullating planes 
at P and Q respectively Q . We approximated 9 by sin 9 where 9 denotes the angle 
between the plane PiPi + iPi + 2 and the plane P,_iPjPj_|_i. If we assign to the plane 
Pj_iPjPj_l_i the equation z — and if we let the segment PjPj+i represent the 
x-axis, then the plane passing through BCD has equation H y — D y z = where D y 
is the projection of the point D on the y-axis. 

Given the equations of two planes Aix + B\y + C\z + D\ =0 and A 2 x + B^y + 
C 2 z + D2 = 0, the angle between the two planes is given by 

A X A 2 + B X B 2 + dC 2 

cos 9 



y/Ai+B'i+CTiy/Ai+B-i + a 

A straightforward computation gives us 
sin6i = ± 



H 2 + D 2 y 



H 

±- 



= ±- 

= ± 
= ± 

(14) = ± 



H 



( e 2 +b 2 — d 2 \ 2 
\ 26 I 



26 

2Hb 



y/4e 2 b 2 - (e 2 + b 2 - d? f 
2Hb 



4A eM 

Hb 
2A e M 
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On the other hand, we can use the local canonical form to approximate the equa- 
tions of the plane P,_jPjPj_|_i and PiPi + iPi + 2 in the Frenet frame coordinates. 
Respectively, we have approximately 

,KT€xS. ,T 



and 



6 ' v 3 

So 



)x + (-)(€!+ e 2 )y + z = 0. 



cos 2 9 

and therefore 



(i + (^( ei -a)( £l + C2 )); 

l + (|( ei -*))a + (I( ei + C2 ))s 



i.e. ^f^- w ±r. Combining this result with Q, we get @. 

With the help of the symbolic computation software Mathematica and using the 
local canonical form of a curve, we computed the Taylor series expansion for both 
expressions. For (O), we obtained 



(15) n =T+^(a-b + 3e) + ^-(b-a + e) + 0(2). 

6k 4 



For (12) we obtained 



(16) T 2 = T+^(a + b + e) + ^-{b-a + e) + 0(2). 

ok 4 



These expansions prove that (|1 1|) and (12) are not equivalent 



There are many ways to approximate r s . For the reasons mentioned above, we 
chose to use (|ll]). In a similar way as for k s , we used ([II]) to obtain the following 
numerically invariant expression for t s . For symmetry reasons and in view of the 
results obtained for n s , we decided on using a centered formula. The result is the 
following five point approximation. 



THP l+1 )-f\P l - 1 ) + {2a + 2b-2d-Zh + g y 1{ l^ ) 
(17) r s (P)-4- 2a + 2b + 2d + h + g 

Where g = d(P,_2, P»-i) as illustrated on figure 1 and k s can be taken to be either 
(|, 1) or (|. 



3.3. Numerical test. We considered the curve a (t) = (cost, sin t, -\/t) whose graph 
is given in figure 12. Note that a is not parametrized by arc length. The signature 
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of this curve is given by the four following quantities. 

_ jofcX a tt \ _ 2Vl6t 3 + 4t 2 + 1 
K " lassl ~ M 3 " (l+4t)3/a ' 

dndt ^ _ 8(8t 2 + 2t-3)Vt 
Ks ~ ~dl~ds ~ if ~ Vl6i 3 + 4t 2 +1(1 + 4t) 3 ' 

_a i 2io«iat« = -2-^(3 + 4t 2 ) 
| a t x a tt | 2 ~~ 16t 3 + 4t 2 + 1 ' 

dr dt _% _ 2(64t 5 - 16i 4 + 240t 3 + 16t 2 - 3) 
Ts ~ ~di ~si ~ if ~~ v / 4TTT(16t 3 + 4t 2 + l) 2 ' 

We tested our numerical expressions on the portion of the curve given by § < t < ^ 
for different partition size. We used k s given by equation (||) for evaluating f s . 
Care was taken to partition the curve very irregularly. In the example presented on 
figures 13-22, we computed the signature for a partition {U}fL built the following 
way 



to 


= 




h 


= t - 


h At 


h 


= h- 


At 

h T 


tz 


= t 2 - 


At 

K T 


U 


= t 3 - 


\- At 


t 5 


= u- 


At 
*" ~2 


t 6 


= t 5 - 


At 
K T etc 



The resulting graphs are presented for At = 0.1, 0.05, 0.025 and 0.0125. For 
simplicity, we graphed two projections of the signature curve. Although a lot of 
the information is lost this way, we believe that it gives a good measure of the 
effectivity of the method. We can see that although the partition is not regular, 
the graphs seem to converge to the exact graphs when At becomes small. Therefore 
we believe that our formulas work properly for generic small partitions. 

4. Conclusion 

We now have in our hands good approximations for the Euclidean and affinc 
signature of a planar curve. The formulas we obtained are invariant under the 
action of the Euclidean and affine group respectively. Another important property 
is that they are valid for any fine partition of the given curve. In a near future, we 
hope to use them for the recognition of planar curves and perhaps also for some 
other applications of the signature curve among the numerous possible ones. In 
particular, we wish to improve the numerical results obtained with noisy images 
and prove that the signature can be a practical tool of object recognition. 

We also have good approximations for the four differential invariants which pa- 
rametrize the Euclidean signature of a space curve, namely the curvature, the tor- 
sion and their derivatives with respect to arc length. The formulas we obtained are 
Euclidean invariants and valid for any fine partition of the given curve, although 
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better results are obtained with equidistant partitions. In particular, we showed in 
section 3 that f = ±6 • -^j^ is a good approximation for the torsion. We believe 
that this expression has a value on its own as it gives an easy visual understand- 
ing of the concept of torsion. Space curve recognition (for example blood vessels 
and trajectory of particles) is an interesting computer vision problem where our 
formulas could find applications. 
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Figure 3. Three Different Partitions of the Initial Curve Used in 
the Planar Euclidean Case. 




a) Regular Partition. 



b) Irregular Partition. 



c) Very Irregular Partition. 



12 



MIREILLE BOUTIN 



Figure 4. Approximations of the Euclidean Signature Curve Ob- 
tained with . 




With Partition 3-a. With Partition 3-b. With Partition 3-c. 



Figure 5. Approximations of the Euclidean Signature Curve Ob- 
tained with 





With Partition 3-a. 



With Partition 3-b. 



With Partition 3-c. 
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FIGURE 6. Approximations of the Euclidean Signature Curve Ob- 
tained with (JiJ). 




With Partition 3-a. With Partition 3-b. With Partition 3-c. 



FIGURE 7. Approximations of the Euclidean Signature Curve Ob- 
tained with (JsJ) . 




With Partition 3-a. With Partition 3-b. With Partition 3-c. 



Figure 8. Approximation of the Euclidean Signature Curve Ob- 
tained with (|gJ) . 




With Partition 3-a. 



With Partition 3-b. 



With Partition 3-c. 
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Figure 9. Exact Data Used for Testing the Planar Affine Case. 




Figure 10. Three Different Partitions of the Initial Curve Used 
in the Planar Affine Case. 




a) Regular Partition. 



b) Irregular Partition. 



c) Very Irregular. 
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FIGURE 11. Approximations of the Affine Signature Curve Ob- 
tained with 
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FIGURE 12. Approximations of the Affine Signature Curve Ob- 
tained with (10). 
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Figure 13. Curve Used to Test the Formulas Proposed for the 
Spatial Signature 




The Curve a(t) = (cost, sint, y/t), for < t < 6tt 



Figure 14. Projections of the Exact Signature for the Curve Above 
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FIGURE 16. Approximations Obtained with At = 0.05 
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Using k and R s ,5 Using k and R s ,4 Using k and k s ,3 
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Figure 17. Approximations Obtained with At = 0.025 




Figure 18. Approximations Obtained with At = 0.0125 
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Figure 19. Approximations Obtained with At = 0.1 
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FIGURE 20. Approximations Obtained with At = 0.05 
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Figure 21. Approximations Obtained with At = 0.025 
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Figure 22. Approximations Obtained with At = 0.0125 
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